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The paper On the multifractal nature of fully developed turbulence and chaotic systems, by R. 
Benzi et al. published in this journal in 1984 ( vol 17, page 3521) has been a starting point of many 
investigations on the different faces of selfsimilarity and intermittency in turbulent phenomena. 
Since then, the multifractal model has become a useful tool for the study of small scale turbulence, 
in particular for detailed predictions of different Eulerian and Lagrangian statistical properties. 
In the occasion of the 50-th birthday of our unforgettable friend and colleague Giovanni Paladin 
(1958-1996), we review here the basic concepts and some applications of the multifractal model for 
turbulence. 



I. INTRODUCTION 



The idea of the multifractal approach to fully developed turbulence has been introduced by Giorgio Parisi and Uriel 
Frisch during the Summer School Turbulence and predictability of geophysical fluid dynamics held in Varenna in June 
1983 One of us (A.V.) had the chance to participate to that school and then to coauthor, with R. Benzi, G. Parisi 
and G. Paladin, the paper published in this journal where the word multifractal appeared for the first time 

From a technical point of view the idea of the multifractal is basically contained in the large deviation theory [3, 0| 
which is an important chapter of the probability theory. However the introduction of the multifractal description in 
the 1980's had an important role in statistical physics, chaos and disordered systems. In particular to clarify in a 
rather net way that the usual idea, coming from the critical phenomena, that just few scaling exponents are relevant, 
is wrong, while an infinite set of exponents is necessary for a complete characterization of the scaling features. 
As pioneering works which anticipated some aspects of the the multifractal approach to turbulence we can cite the 
lognormal theory of Kolmogorov 0] , the contributions of Novikov and Stewart Q and Mandelbrot ■ 

This paper has no pretense to be a survey of the many applications of the multifractal description in chaos, disordered 
systems and natural phenomena; for general reviews on these aspects see [1, • For a more mathematically oriented 
treatment see [ill]. Our aim is a discussion on the the use of the multifractal methods in the study of the scaling 
features of fully developed turbulence. 

The paper is organized as follows. Section |TI] is devoted to the introduction of the multifractal model of turbulence 
and its connections with the f{a) vs a formalism introduced by Halsey et al. 12], and the large deviations theory [1,0]- 
In Section IIIII we discuss the implication of multifractality on Eulerian features, namely the statistical properties of 
the velocity gradients and the existence of an intermediate dissipative range. Section HVl is devoted to the implications 
of the multifractal nature of turbulence on Lagrangian statistics. In Section 5 we present the Lagrangian acceleration 
statistics. Section 6 treats the relative dispersion, in particular we discuss the multifractal generalization of the 
classical Richardson theory. Section 7 is devoted to the multifractal analysis of the dispersion in two-dimensional 
convection. 



II. FROM KOLMOGOROV TO MULTIFRACTALS 



Let us consider the Navier-Stokes equations for an incompressible fluid: 

dv 1 

— -|-(v-V)v = — Vp + i^Av + F, V-v = 0. (1) 
ot p 

Because of the nonlinear structure of the equation, an analytical treatment is a formidable task. For instance in the 
3D case a theorem for the existence of global solution for arbitrary v is still missing. 

For a perfect fluid (i.e. = 0) and in absence of external forces (F = 0), the evolution of the velocity fleld is given by 
the Euler equation, which conserves the kinetic energy for smooth solutions. In such a case, introducing an ultraviolet 
cutoff Kmax On the wave numbers, it is possible to build up an equilibrium statistical mechanics simply following 
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the standard approach used for the Hamiltonian statistical mechanics. On the other hand, because of the so called 
dissipative anomaly p^fl^. in 3D the limit — > is singular and cannot be interchanged with Kj^ax ~* oOj therefore 
the statistical mechanics of an inviscid fluid has a rather limited relevance for the Navier-Stokes equations at very 
high Reynolds numbers {Re — VL/iy, where V and L are the typical speed and length of the system, respectively). 

In addition, mainly as a consequence of the non-Gaussian statistics, even a systematic statistical approach, e.g. in 
term of closure approximations, is very difficult [isL [l^. 

In the fully developed turbulence (FDT) limit, i.e. v ^ 0, and in the presence of forcing at large scale, one has a 
non equilibrium statistical steady state, with an inertial range of scales, where neither energy pumping nor dissipation 
acts, which shows strong departures from the equipartition (isl . [Ti| . A simple and elegant explanation of the main 
statistical features of FDT is due to Kolmogorov [l^ : in a nutshell, it is assumed the existence of a range of scales 
where the energy - injected at the scale L - flows down (with a cascade process, as remarked by Richardson [l^) to 
the dissipative scale £d ~ LRe~^^'^, where it is dissipated by molecular viscosity. Since, practically, neither injection 
nor dissipation takes place in the inertial range, the only relevant quantity is the average energy transfer rate e. 
Dimensional counting imposes the power law dependence for the second order structure function (SF) 

S2{1) - {dv}) = {{V{X V{x)f) OC £-2/3^2/3 ^ (2) 

where, for the sake of simplicity, we ignore the vectorial nature of the velocity field. The scaling law ([2|) is equivalent to 
a power spectrum E{k) oc e^^^k^^^^ in good agreement with the experimental observations. The original Kolmogorov 
theory (often indicated as K41) assumes self-similarity of the turbulent flow. As a consequence, the scaling behavior 
of higher order structure functions 

Sp{£)^{\v{x + £)-v{x)\P)^6'' (3) 
is described by a single scaling exponent: C,p — p/3. 

A. The multifractal model 

The Navier-Stokes equations arc formally invariant under the scaling transformation: 

with A > (indeed the Reynolds number VL/v is invariant under the above transformations). 

The exponent h cannot be determined with only symmetry considerations, nevertheless there is a rather natural 
candidate: h — 1/3. Such a value of the exponent is suggested by the dimensional argument of the K41 and also, 
and more rigorously, by the so-called "4/5 law", an exact relation derived by Kolmogorov from the Navier-Stokes 
equations [lH, , which, under the assumption of stationarity, homogeneity and isotropy, states 

{Svf\{e))^~^e£, (4) 

where Sv^{£) is the longitudinal velocity difference between two points at distance £. 

We can say that the K41 theory corresponds to a global invariance with h = 1/3 and therefore (p = p/3. This 
result is in disagreement with several experimental investigations ^13i . il8l] which have shown deviations of the scaling 
exponents from p/3. This phenomenon, which goes under the name of intermittency [l3| . is a consequence of the 
breakdown of self-similarity and implies that the scaling exponents cannot be determined on a simple dimensional 
basis. 

A simple way to modify the K41 consists in assuming that the energy dissipation is uniformly distributed on 
homogeneous fractal with dimension Dp < 3. This implies 6ve{x) ~ i'^ with h = {Dp — 2)/3 ior x on the fractal and 
dve{x) non singular otherwise. This assumption (called absolute curdling or /3- model) gives: 

(p^^^P+{5-Dp). (5) 

Such a prediction, with Dp ^ 2.83 , is in fair agreement with the experimental data for small values of p, but higher 
order scaling exponents give a clear indication of a non linear behavior in p. 

The multifractal model of turbulence [l|, [l^ [l^ assumes that the velocity has a local scale-invariance, i.e. there is 
not a unique scaling exponent h such that 5vg ~ , but a continuous spectrum of exponents, each of which belonging 
to a given fractal set. In other words, in the inertial range one has 



(6) 
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if X S Sh, and Sh is a fractal set with dimension D{h) and h G {hmim hmax)- The probabihty to observe a given 
scahng exponent h at the scale i is Peih) ~ £3-D(h) therefore one has 



For £ 1, a steepest descent estimation gives 

Cp = mm{hp + 3- D{h)} = h*p + 3- D{h*) 



(7) 



(8) 



where h* = h*{p) is the solution of the equation D'{h*{p)) — p. The Kolmogorov "4/5" law dU imposes Ca = 1 which 
implies that 



D{h) < 3/1 + 2, 



(9) 



with the equality realized by /i*(3). The Kolmogorov similarity theory corresponds to the case of only one singularity 
exponent h — 1/3 with D{h — 1/3) — 3. 

Of course the computation of D{h), or equivalently C,p, from the NSE is not at present an attainable goal. A first 
step is a phenomenological approach using multiplicative processes. Let us briefly remind the so called random /3- 
model t^. This model describes the energy cascade in real space looking at eddies of size in = 2~"L, with L the 
length at which the energy is injected. At the n- th step of the cascade a mother eddy of size in splits into daughter 
eddies of size in+i, and the daughter eddies cover a fraction /3j (0 < f3j < 1) of the mother volume. As a consequence 
of the fact that the energy transfer is constant throughout the cascade one has for the velocity differences w„ = Sv^^ 
on scale in is non negligible only on a fraction of volume Ylj Pj , and it is given by 



= voil/'l[(^j 



-1/3 



(10) 



where the /3j's are independent, identically distributed random variables. Phenomenological arguments suggest: Pj = 1 
with probability x and (3j = B = 2~(^~^'''"'") with probability 1 — a; The above multiplicative process generates a 



two-scale Cantor set, which is a rather common structure in chaotic systems. The scaling exponents are: 



corresponding to 



D{h) = 3 + (3/i- l) 



l + ln- 



/I - 3/i\ 
\l-x) 



ih Ino 



3/i 



(11) 



(12) 



The two limit cases are x = 1, i.e. the K41, and x = which is the f3- model with Dp = 2 + 3/imm- Using x = 7/8, 
hmin = (i.e. B = 1/2) one has a good fit for the (p of the experimental data at high Reynolds numbers. 

Of course it is not so astonishing to find a model to fit the experimental data. Indeed, there are now many 
phenomenological models for D{h) which provide scaling exponents in agreement with experimental data. A popular 
one is the so-called She-Leveque model [20J which is reproduced by the multifractal model with 



D{h) = 1 



2/3 -3/1-1 



1 - In 



/2/3 - 1 - 3/i 
V 21n/3 



and gives for the scaling exponents 



Cp 



2/3-1 



p + 2(1-/3^/3) 



(13) 



which are close to the experimental data for /3 = 2/3. Another important model, which was introduced by Kolmogorov 
himself without reference to the multifractal model, is the log-normal model which will be discussed in Section fll Bl 
The relevance, and the success, of the multifractal approach is in the possibility to predict, and test, nontrivial 
statistical features, e.g. the pdf of the velocity gradient, the existence of an intermediate dissipative range and precise 
scaling for Lagrangian quantities. Once D{h) is obtained by a fit of the experimental data, e.g. from the ^p, then all 
the predictions obtained in the multifractal model framework must be verified without additional free parameters. 
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FIG. 1: Structure function scaling exponents (p plotted vs p. Circles and triangles correspond to the data of Anselmet at al. 
[l^ . The solid line corresponds to Kolmogorov scaling p/3; the dashed line is the random beta model prediction with 
B — 1/2 and x — 7/8; the dotted line is the She-Leveque prediction (|13p with P = 2/3. 



B. Relation between the original multifractal model and the /(a) vs a 



The multifractal model for the FDT previously discussed is linked to the so called f{a) vs a description of the 
singular measures (e.g. in chaotic attractors) [l3,|Tl,|ll|. In order to show this connection let us recall the Kolmogorov 
revised theory Q (called K62) stating that the velocity increments 5vi scales as {ei£y^^ where ee is the energy 
dissipation space-averaged over a cube of edge i. Let us introduce the measure /i(x) = e(x)/e, a partition of non 
overlapping cells of size £ and the coarse graining probability 



where A/(xi) is a cube of edge i centered in x^, of course eg ~ i^^P{£). Denoting with a the scaling exponent of Pi 
and with /(a) the fractal dimension of the subfractal with scaling exponent a, we can introduce the Renyi dimensions 
dr,: 



where the sum is over the non empty boxes. A simple computation gives 

{p ~ l)dp = min[pa - f{a)] . 

a 

Noting that < >= ^■^ ^ we have 

< £(P- -3) 

therefore one has the correspondence 

a -2 



D{h)^f{a) , Cp-^ + (f -l)(dz-3) 



Of course the result C3 = 1, once assumed 5v ^ [eit)^/^, holds for any /(a). Let us note that the lognormal theory 
K62 where 
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is a special case of the multifractal model, where there are not restriction on the values of h and D{h) is a parabola 
with a maximum at Dp = 3: 

D{h) = -—h^ + -{2 + fi)h 



2n 2 ' " 8/Lt 
and the parameter jj, is determined by the fluctuation of Ine^. 

C. A technical remark on multifractality 

To obtain the scaling behavior of Sp{t) ~ {1/ L)'^^ given by ([7]) with obtained from ([5]), one has to assume that the 
exponent p/i + 3 — D{h) has a minimum, (p, which is a function of /i, and that such an exponent behaves quadratically 
with h in the vicinity of the minimum. This is the basic assumption to apply the Laplace's method of steepest descent 
[2^ . The point wc would like to recall here is that, for small separations, £, it is true that Sp{() ~ (£/L)^p but with 
a logarithmic prcfactor: 



-1/2 



(14) 



Such a prefactor is usually not considered in the naive application of Laplace method leading to (O . The presence of 
such logarithmic correction, if present, would clearly invalidate the 4/5-th law ([3]), one of the very few exact results 
in fully developed turbulence. 

The question on whether such logarithmic correction is likely has quantitatively been addressed by Frisch et al. [23j . 
There, exploiting the refined large-deviations theory, the Authors were able to explain in which way the logarithmic 
contribution cancels out thus giving rise to a prediction fully compatible with the naive (a priory unjustified) procedure 
to extract the scaling behavior (O. The key point is that the leading order large deviation result for the probability 
Pi{h) to be within a distance i of the set carrying singularities of scaling exponent between h and h -f dh, 

3~D{h) 



Piih) ^[j] , (15) 



must be generalized to take into account next subleading order. In doing so, as a result one obtains [22 







1/2 


^ 3-D(h) 


[-mil 






L 





(16) 



which contains subleading logarithmic correction. It is worth observing that despite the multiplicative character of the 
logarithmic correction one speaks of "subleading correction" . This is justified by the fact that the correct statement 
of the large-deviations leading-order result involves the logarithm of the probability divided by the logarithm of the 
scale. The correction is then a subleading additive term. 
Once the expression (jl6p is plugged in the integral 

Sp{e)^ J dhp,{h) (^ly (17) 

and the saddle point estimation is carried out according to [l^, logarithms disappear and the expected 4/5-th law 
emerges. 

It is worth mentioning that the presence of a square root of a logarithm correction in the multifractal probability 
density had already been proposed by [2^ on the basis of a normalization requirement. In that paper, the Authors 
observed that without such a correction the singularity spectrum /(a) comes out wrong; they also pointed out that 
a similar correction has been proposed by [2^ in connection with the measurement of generalized Renyi dimensions. 

III. IMPLICATIONS OF MULTIFRACTALITY ON EULERIAN FEATURES 

At first we note that a consequence of presence of the intermittency the Kolmogorov scale does not take a unique 
value. The local dissipative scale is determined by imposing the effective Reynolds number to be of order unity: 



(18) 



6 



therefore the dependence of Ijj on ft, is thus 

£D{h) - LRe-^ (19) 

where Re = Re{L) is the large scale Reynolds number (26j . 

In this section we will show that the fluctuations of the dissipative scale, due to the intermittency in the turbulent 
cascade, is relevant of the statistical features of the velocity differences and velocity gradient, and in addition it implies 
the existence of an intermediate region between the inertial and dissipative range [27| . 



A. The Pdf of the velocity differences and velocity gradient 

Let us denote by s the longitudinal velocity gradient. Such a quantity can immediately be expressed in terms of 
the singularity exponents h as 



Svep ffh-i 



i;o^^-i = «o'^'V^ (20) 



where we used the fact that Sv^ ^ vqI^ from ^ and we have exploited p^ . From ([^0]) we realize that we can 
easily express the probability density function (PDF) of s (for a fixed ft), Ph{s), in terms of the PDF, n(Vo), of the 
large-scale velocity differences Vq, with vq = |Vb|. The latter PDF is indeed known to be accurately described by the 
Gaussian distribution [111. The link between the two PDFs is given by the standard relation: 



p,.(s) = n(Fo) 

from which one immediately gets: 



ds 



(21) 



V 2 TTTTTT 



Ph{s)r^— e ^<V> . (22) 



s 



The K41 theory and the f3 model correspond to ft = 1/3 and ft = {Dp — 2)/3, respectively. In both cases, a 
stretched exponential form for the PDF is predicted with an exponent, 1 + ft, larger than one. Experimental data (see 
e.g. j2^ |30|) are not consistent with such a prediction being actually compatible with a stretched exponent whose 
value is smaller than one. 

The multifractal description has thus to be exploited to capture those experimental evidences. To do that we recall 
the expression (fTO]) for the random /? model 

n 

Vn = voel/'^l[f3-'/' (23) 
from which the probability distribution of the velocity increments Vn reads: 

p r n n 

PM = / U{Vo)dVo I S{vn - voll'^ n /^r'^') n PMW^- (24) 

1=1 j=i 

Here niPi) is the probability density of the PiS assumed to be of the form: 

^i{^3^) = x(5(A - 1) + (1 - x)5{l3, - B) (25) 

with B = 2-(i-3''"""). 

Since the /3j's are identically distributed, the above integral becomes: 



PM=j2 (")a:""^(l-x)^B4^/3^-i/3g-^^'"'''^''"''" 



(26) 



with C = (2(Fq^}) ^. It is easy to see [3l| the passage of the above PDF from a Gaussian form at large scales (small 
n) to an exponential- like form at small scales (large n). 
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To obtain the gradient PDF from pH)) it is sufficient to stop the sum at n = iV such that vnIn/v = 1. This is 
equivalent to say 



1% ^ 2 



-2N 



or N = 



In- 



s 21n2' 

By noting that B^^ — (i/ / s)^^^'^""" the resuhing gradient PDF reads: 



N 



(l+29)/3 



-Cv 



(2 + <j)/3|„|(4~,)/3 



(27) 



(28) 



where q = K{1 — 3hmin)/N. The K41 prediction corresponds to considering only the term K = with x = 1. 
We already discussed that x = 7/8 and hmin = provides a good fit for the scaling exponents Cp of the structure 
functions in the limit of high Reynolds numbers. The same parameters give a PDF behavior in good agreement with 
available experimental data (see [3l| and Figure [2]). 
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FIG. 2: Log-linear plot of the PDF of velocity gradients s rescaled with the rms value. Points represent experimental data 
from ^30j, solid line is the multifractal prediction with the random beta model, dotted and dashed lines represent the K41 and 
beta model results respectively. 



B. Intermediate dissipative range 



Let us now show, that as consequence of the fluctuations of the dissipative scale one has the existence of an 
intermediate region (the Intermediate Dissipative Range, IDR) between the inertial and dissipative range [13]. The 
presence of fluctuations of £d, see (IT^ . modifies the evaluation of the structure functions ([7]): for a given £, the saddle 
point evaluation remains unchanged if, for the selected exponent h*{p), one has £D{h*{p)) < L If, on the contrary, 
the selected exponent is such that £D{h*{p)) > I the saddle point evaluation is not consistent, because at scale I 
the power-law scaling ([6]) is no longer valid. In this intermediate dissipation range the integral in ([7]) is dominated 
by the smallest acceptable scaling exponent h(C) given by inverting (|19p . and the structure function of order p a 
pseudo-algebraic behavior, i.e. a power law with exponent ph[C) + 3 — D{h{£)) which depends on the scale £. Taking 
into account the fluctuations of the dissipative range [27] , one has for the structure functions 

I gHi)p+3-Dihie)) if i^(hrmn) < £ < ^d(/i*(p)) ■ ^ ' 
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A simple calculation [13|, |27| shows that it is possible to find a universal description valid both in the inertial and in 
the intermediate dissipative ranges. Let us discuss this point for the energy spectrum E{k). Introducing the rescaled 
variables 

InEik) , „ Ink , , 

^(^) = p and e = —— (30) 



In Re In Re 



one obtains the following behavior 




nO)-^ o on , eD{e^^-l) for j^<e<j^ ^^'> 

The prediction of the multifractal model is that ln£'(fc)/lni?e is an universal function of lnfc/lni?e. This is in 
contrast with the usual scaling hypothesis according which \nE{k) should be a universal function of \n{k/k]:))). The 
multifractal universality has been tested by collapsing energy spectra obtained from turbulent flow in a wide range of 
Re fal, see also (sl. 



C. Exit times for turbulent signals and the IDR 



In the following we will discuss a method alternative to the study of the structure functions which allows for a 
deeper understanding of the IDR. 

Basically in typical experiments one is forced to analyze one-dimensional string of data v{t), e.g. the output 
of hot-wire anemometer, and the Taylor Frozen- Turbulence Hypothesis is used to bridge measurements in space 
with measurements in time. As a function of time increment, r, structure functions assume the form: Spir) =< 
[{v{t + t) — v{t)]^>. In the inertial range, td <C r <^ Tq (where Tq = Lq/Vo, and the dissipative time, td = £d/Vo) 
the structure functions develop an anomalous scahng behavior: S'p(t) ~ r'»p, where r ~ i/Vo. 

The main idea, which can be applied both to experimental and synthetic data, is to take a time sequence v(t), and 
to analyze the statistical properties of the exit times from a set of defined velocity-thresholds. More precisely, given 
a reference initial time to with velocity v{to), we define t{Sv) as the first time necessary to have an absolute variation 
equal to Sv in the velocity data, i.e. \v{to) — v{to + t{6v))\ = Sv. By scanning the whole time series we recover the 
probability density functions of T{dv) at varying 6v from the typical large scale values down to the smallest dissipative 
values. Positive moments of T{dv) are dominated by events with a smooth velocity field, i.e. laminar bursts in the 
turbulent cascade. Let us define the Inverse Structure Functions (Inverse-SF) as [33l. [33|: 

Sp((5w) =<tP{6v)> . (32) 

It is necessary to perform weighted average over the time-statistics in a weighted way. This is due to the fact that by 
looking at the exit-time statistics we are not sampling the time-series uniformly, i.e. the higher the value of t{Sv) is, 
the longer it remains detectable in the time series. 



It is possible to show [35| that the sequential time average of any observable, A, based on exit-time statistics, (^)e, 



is connected to the uniformly-in-time multifractal average by the relation: 



For A — tP{6v) the above relations becomes: 



(^) = ^ . (33) 

{T)e 



According to the multifractal description we assume that, for velocity thresholds corresponding to inertial range values 
of the velocity differences the following dimensional relation is valid: 

S-,v ^ t{5v) - 5v^/'' , 

and the probability to observe a value r for the exit time is given by inverting the multifractal probability, i.e. 
P(t ~ Sv'^/^) ~ (57;[3-^('»)]/'»_ With this ansatz in the inertial range one has: 

T.p{5v) - / dh 5v[p+3--D('')l/'' ^ JvXp (35) 
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where with the Laplace method one obtains: 

Xp = min{b + 3-i?(M]/M • (36) 

h 

Let us now consider the IDR properties. 

For each p, the saddle point evaluation selects a particular h — hs{p) where the minimum is reached. Let us also 
remark that from ()35|1 we have an estimate for the minimum value assumed by the velocity in the inertial range given 
a certain singularity h: Vmih) ~ (5^^(^)W ~ i/h/(i+h)^ Therefore, the smallest velocity value at which the scaling ([55]) 
still holds depends on both i> and h. Namely, 6vmip) i/''^ (p). The most important consequence is that for 
Sv < 5vm{p) the integral (|35p is not any more dominated by the saddle point value but by the maximum h value 
still dynamically alive at that velocity difference, l/h{5v) = — 1 — \og{v) / \og{5v) . This leads for 5v < Svjn{p) to a 
pseudo-algebraic law: 

j:,isv)r.svip+^-D(''(^-))y''(s-y (37) 

The presence of this p-dependent velocity range, intermediate between the inertial range, Sp((5w) ~ Sv^p, and the 
dissipative scaling, T,p{6v) ~ 6v^, is the IDR signature. 




FIG. 3: Inverse Structure Functions 'Ei{5v). The straight lines shows the dissipative range behavior (dashed) Ei(5u) --^ 5v, and 
the inertial range non intermittent behavior (dotted) Tji{Sv) ~ (Sv)^. The inset shows the direct structure function Si(t) with 
superimposed the intermittent slope (i — 0.39. 

In Figure [3] we show Y.i{Sv) evaluated on a string of high- Reynolds number experimental data as a function of the 
available range of velocity thresholds Sv. This data set has been measured in a wind tunnel at Re\ ^ 2000. One can 
see that the scaling is very poor. On the other hand, (inset of Figure [3]), the scaling behavior of the direct structure 
functions < |(5w(t)| >^ t'-i is quite clear in a wide range of scales. This is a clear evidence of IDR's contamination 
into the whole range of available velocity values for the Inverse-SF cases. 

Let us now go back to the statistical properties of the IDR. In order to study this question we have smoothed the 
stochastic synthetic field, v{t) (see Appendix) by performing a running-time average over a time- window, ST. Then 
we compare Inverse-SF obtained for different Reynolds numbers, i.e. for different dissipative cut-off: Re ^ ST~^/^. 
The expression (|57)) predicts the possibility to obtain a data collapse of all curves with different Reynolds numbers 
by rescaling the Inverse-SF as follows [13, [36| : 

- ln(Sp((5w))/ln((5T/5To) vs. -\n{5v/U)/\n{ST/STo) , (38) 

where U and 8Tq are adjustable dimensional parameters. 

Figure |4] shows the rescaling (|38p of the Inverse-SF, Ei((5w), both for the synthetic field at different Reynolds 
numbers and for the experimental signals. As it is possible to see, the data-collapse is very good. This is a clear 
evidence that the poor scaling range observed in Figure 0] for the experimental signal can be explained as the signature 
of the IDR. 
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FIG. 4: Data collapse of the Inverse-SF, Ei((5«), obtained by the rescaling (|38|l for the smoothed synthetic signals (with time 
windows: 5T — 4.8 • 10^'' , 3 • 10~^ , 2 ■ 10"") and the experimental data (EXP). The two straight lines have the dissipative 
(solid line) and the inertial range (dashed) slope. 



IV. THE RELATION BETWEEN EULERIAN AND LAGRANGIAN STATISTICS 

A problem of great interest concerns the study of the spatial and temporal structure of the so-called passive fields, 
indicating by this term quantities transported by the flow without affecting the velocity field. The paradigmatic 
equation for the evolution of a passive scalar field 9{x, t) advected by a velocity field v(x, t) is [s^l 

946* + V ■ (v 6*) = X V^6' (39) 

where x is the molecular diffusion coefficient. 

The problem can be studied through two equivalent approaches, both due to Euler [s^. The first, referred to 
as "Eulerian" , deals at any time with the field 9 in the space domain covered by the fluid; the second considers the 
time evolution of trajectories of each fluid particle and is called "Lagrangian" . 

The motion of a fluid particle is determined by the differential equation 

^=v(x,t) (40) 

which also describes the motion of test particles, for example a powder embedded in the fluid, provided that the 
particles are neutral and small enough not to perturb the velocity field, although large enough not to perform a 
Brownian motion. Particles of this type are commonly used for fiow visualization in fiuid mechanics experiments [3^ . 
We remark that the complete equation for the motion of a material particle in a fluid when density and volume effects 
are taken into account can be rather complicated [13, |4lj . 

The Lagrangian equation of motion (|40p formally represents a dynamical system in the phase space of physical 
coordinates. By very general considerations, it is now well established that even in regular velocity field the motion 
of fiuid particles can be very irregular (4^.l43j. In this case initially nearby trajectories diverge exponentially and one 
speaks of Lagrangian chaos or chaotic advection. In general, chaotic behaviors can arise in two-dimensional flow only 
for time dependent velocity fields, while it can be present even for stationary velocity fields in three dimensions. 

If X = 0, it is easy to realize that ([M)) is equivalent to (^0]) . Indeed, we can write 

6i(x,<) = 6i(T-*x,0) (41) 

where T is the formal evolution operator of (|40|) : x(t) ~ r*x(0). 

Taking into account the molecular diffusion x, ([M]) is the Fokker-Planck equation of the Langevin equation [il] 

^=v(x,i)+r;(i) (42) 
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where ry is a Gaussian process with zero mean and variance 

{v^(t)7^,it')) = 2xS^JS{t-t'). (43) 



The dynamical system (|40|) becomes conservative in the phase space in the case of an incompressible velocity field 
for which 

V • V = (44) 

In two dimensions, x — (xi,X2), the constraint ()44[) is automatically satisfied by introducing the stream function 
V'(x,t) 

^V- dip 
vi = ^ — , V2 = - — (45) 
0x2 0x1 



and the evolution equation becomes 



dxi _ dtp dx2 _ dtp 

dt dx2 ' dt dxi ' 



i.e. formally a Hamiltonian system with the Hamiltonian given by the stream function ■0. 

The presence of Lagrangian chaos in regular flows is a remarkable example of the fact that, in general, it is very 
difficult to relate Lagrangian and Eulerian statistics. For example, from very complicated trajectories of buoys one 
cannot infer the time-dependent circulation of the sea. In the following Sections we will see that in the case of fully 
developed turbulence, the disordered nature of the flow makes this connection partially possible at a statistical level. 

The equation of motion (PU]) shows that the trajectory of a single particle is not Galilean invariant, i.e. invariant 
with respect to the addition of a mean velocity. The most general Galilean invariant statistics, which is ruled by 
small scale velocity fluctuation, is given by multi-particle, multi-time correlations for which we could expect universal 
features. In the following we will consider separately the two most studied statistics: single-particle two-time velocity 
differences and two-particle single-time relative dispersion. 



A. Single particle statistics: multifractal description of Lagrangian velocity differences 

The simplest Galilean invariant Lagrangian quantity is the single particle velocity increment 5w{t) — v(fo +0 ~'^(^o) j 
where v(i) = v(x(t),t) denotes the Lagrangian velocity of the particle at x(t) and the independence on to is a 
consequence of the stationarity of the flow. Dimensional analysis in fully developed turbulence predicts p^ . |45| 

{5v,{t)6vj{t)) = CoetS^, (47) 

where e is the mean energy dissipation and Co is a numerical constant. The remarkable coincidence that the variance 
of 5v(t) grows linearly with time is the physical basis on which stochastic models of particle dispersion are based. It 
is important to recall that the "diffusive" nature of (|47p is purely incidental: it is a direct consequence of Kolmogorov 
scaling in the inertial range of turbulence and is not directly related to a diffusive process (i.e. there is no decorrelation 
justifying the applicability of central limit theorem). 

Let us recall briefly the argument leading to the scaling in ([TT]) . We can think at the velocity v{t) advecting the 
Lagrangian trajectory as the superposition of the different velocity contributions coming from turbulent eddies (which 
also move with the same velocity of the Lagrangian trajectory). After a time t the components associated to the 
smaller (and faster) eddies, below a certain scale € are decorrelated and thus at the leading order one has 5v{t) ~ 6v{(!.). 
Within Kolmogorov scaling, the velocity fluctuation at scale £ is given by 5v{(P) ^ Vq{1/ L)^/'^ where Vb represents the 
typical velocity at the largest scale L. The correlation time of 5v{t) scales as t{1) ~ tq{E/ LY^^ and thus one obtains 
the scahng in (|T7|) with e = Vq /tq. 

Equation (|47p can be generalized to higher order moments with the introduction of a set of temporal scaling 
exponents 

{5v{tY) ^ (t/ro)«'' (48) 

The dimensional estimation sketched above gives the prediction = p/2 but one may expect corrections to the 
dimensional scaling in the presence of intermittency. 
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A generalization of the above results which takes into account intermittency corrections can be easily developed by 
using the multifractal model (46l. |47|. The dimensional argument is repeated for the local scaling exponent h, giving 
dv{t) - Vo{t/To)''/'-^-''\ Integrating over the h distribution one ends with 

[ph-D{h)+3]/{l-h) 



{Svm 



^0 



dh ( - 

To 



(49) 



where D{h) is the Eulerian fractal dimension (i.e. 
(■(g) = mmh[qh — D{h) + 3]). In the limit t/ro 
giving the prediction 



related to the Eulerian structure function scaling exponents by 
0, the integral can be estimated by a steepest descent argument 



mm 

h 



ph-D{h)+3 



1 - h 



(50) 



The standard inequality in the multifractal model D{h) < 3h + 2 implies for ((50l) that even in the presence of 
intermittency ^2 = 1- Physically, this is a consequence of the fact that energy dissipation is raised to the first power, 
in gZl). ^ 

Experimental results 48] have shown that even at large Reynolds number the scaling (|¥71) is not clearly observed. 
Therefore the dimensionless constant Cq is known with large uncertainty, if compared with the Kolmogorov constant. 

Intermittency in Lagrangian velocity differences is evident by looking at the pdf of 5v{t) at different time lags, as 
shown in Fig. [S] For large time delays the pdf are close to Gaussian while decreasing t they develop larger and larger 
tails, implying the breakdown of self-similarity. 




FIG. 5: Probability density functions of velocity increments for a numerical simulation at R\ — 284. Curves refer to time 



increments t = (97, 24, 6, 0.7)rD from inside to outside (and to the acceleration), 
function of time delay is shown. Figure from 49] . 



In the inset kurtosis {5v{t)^) / (6v{ty}^ 



Higher order Lagrangian structure functions are shown in Fig.[6]for a set of direct numerical simulations at Rx = 284 
[jgl . [50 |. Despite the apparent scaling observed in the log-log plot, the computation of local slopes does not give a 
definite value of scaling exponents. Assuming ^2 = 1 as predicted by (1501) . one can measure the relative scaling 
exponent ^p/^2 by using the so-called extended self-similarity procedure [51|. As shown by the inset of Fig. [6l we 
observe a well defined scaling in the range of separations lOr^i < t < SOtjj- The values of the relative exponents 
estimated with this method, £,4/^,2 = 1.7 ± 0.05, 1^5/^2 = 2.0 ± 0.05, £,6/£,2 = 2.2 ± 0.07, are in good agreement with 
those predicted by the multifractal model ([50]) . 

The multifractal prediction ((50|) has been checked also in simplified Lagrangian model based on the shell model of 
turbulence 14711. 



V. LAGRANGIAN ACCELERATION STATISTICS 



Acceleration in fully developed turbulence is an extremely intermittent quantity which display fluctuations up to 
80 times its root mean square [53 |. These extreme events generate very large tails in the pdf of acceleration which 
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FIG. 6: Log-log plot of Lagrangian structure functions of order 2, 4 and 6 (bottom to top) for the three components of 
the velocity obtained from DNS at R\ = 284. Small differences between the components are observed at large delays as a 
consequence of small anisotropies at large scales. In the insets the local slopes (bottom right) and relative slopes (top left) 
are shown for the x-component. Relative slopes are defined as dSp{t)/dS2{t) for p = 4,6. The horizontal lines represent the 
multifractal prediction. Figure from [49l |. 



are therefore expected to be very far from Gaussian. 

We remark that even within non-intermittent Kolmogorov scaling, acceleration pdf is expected to be non-Gaussian. 
Indeed acceleration can be estimated from velocity fluctuations at the Kolmogorov scale as 

a = ^ (51) 

where td = r]/Sv{r]) and the Kolmogorov scale 77 is given by the condition •q5v['q)/v = 1. By assuming the scaling 
6v{£) ~ Vo{£/L)^ (with h — 1/3 for Kolmogorov scaling) one obtains 



and therefore 



1 (Y^ 



l-2h 
1 + h 



(52) 



(53) 



Assuming a Gaussian distribution for large scale velocity fluctuations Vq (which is, as already observed, consistent 
with many experimental and numerical observations), and taking h = 1/3, one obtains for the pdf of a a stretched 
exponential tail p{a) ^ exp(— Ca®/^). 

In the presence of intermittency the above argument has to be modifled by taking into account the fluctuations 
of scaling exponent. In the recent years, several models have been proposed for describing turbulent acceleration 
statistics, on the basis of different physical ingredients. In the following we want to show that the multifractal model 
of turbulence, when extended to describe fluctuation at the dissipative scale, is able to predict the pdf of acceleration 
observed in simulations and experiments with high accuracy [53|. Moreover, as in the case of Lagrangian structure 
functions, the model does not require the introduction of new parameters, a part the set of Eulerian scaling exponents. 
In this sense, multifractal model become a predictive model for Lagrangian statistics. 

The introduction of intermittency in the above argument is simply obtained by weighting ([55]) with both the 
distribution of Vq (still assumed Gaussian, as intermittency is not expected to affect large scale statistics) and the 
distribution of scaling exponent h which can be rewritten as 

7/\3--D('i) fVoL 



Pih) - i f) -{ — ] (54) 
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The final prediction, when written for the dimensionless acceleration a = a/ (a^)^/^, becomes [s^ 

p{a) ^ ^ a[h-^+Dm/^RyW exp (^_ia2(i+,.)/3^^(M^ (55) 

where y{K) ^ x{h - ^ + D{h))/Q + 2{2D{h) + 2h ~ 7)/3 and z{h) = x(l + /i)/3 + 4(2/i - l)/3. The coefHcient 
X is the scaling exponent for the Reynolds dependence of the acceleration variance, (a^) ^ R^, given by % = 
sup^ {2{D(h) — Ah — 1)/(1 + h)). For the non-intermittent Kolmogorov scaling {h = 1/3 and D(l/3) — 3) one obtains 
X — ^ E^nd (|55p recovers the stretched exponential prediction discussed above. 

We note that (j55p may show an unphysical divergence for a ^ for many multifractal models of D(h) at small h. 
This is not a real problem for two reasons. First, the multifractal formalism cannot be extended to very small velocity 
and acceleration increments because it is based on arguments valid only to within a constant of order one. Thus, 
it is not suited for predicting precise functional forms for the core of the pdf. Second, small values of h correspond 
to very intense velocity fluctuations which have never been accurately tested in experiments or by DNS. The precise 
functional form of D(h) for those values of h is therefore unknown. 

In Fig. ([7]) we compare the acceleration pdf computed from the DNS data at R\ — 280 with the multifrac- 
tal prediction ([55]) using for D{h) an empirical model which fits well the Eulerian scaling exponents [20[. The 
large number of Lagrangian particles used in the DNS (see [H^ for details) allows us to detect events up to 80 Ca- 
The accuracy of the statistics is improved by averaging over the total duration of the simulation and all directions 
since the flow is stationary and isotropic at small scales. Also shown in Fig. ([7]) is the non- intermittent prediction 
exp (— As is evident from the figure, the multifractal prediction captures the shape of the 
acceleration pdf much better than the K41 prediction. What is remarkable is that ([55)1 agrees with the DNS data well 
into the tails of the distribution - from the order of one standard deviation a a up to order TOua. We emphasize that 
the only free parameter in the multifractal formulation of p{a) is the minimum value of the acceleration, amin, here 
taken to be 1.5. In the inset of Fig. ([7]) we make a more stringent test of the multifractal prediction ([55)1 by plotting 
a^pia) and which is seen to agree well with the DNS data. 




10 20 30 40 50 60 70 80 90 



FIG. 7: Log-linear plot the acceleration pdf. Points are the DNS data, the solid line is the multifractal prediction and the 
dashed line is the K41 prediction. The DNS statistics were calculated along the trajectories of two million particles amounting 
to 3.6 X 10^ events in total. The statistical uncertainty in the pdf was quantified by assuming that fluctuations grow like the 
square root of the number of events. Inset: a*p(a) for the DNS data (crosses) and the multifractal prediction. 



VI. RELATIVE DISPERSION IN TURBULENCE 

Relative dispersion of two particles is historically the first issue quantitatively addressed in the study of fully 
developed turbulence. This was done by Richardson, in a pioneering work on the properties of dispersion in the 
atmosphere in 1926 [s^l, and then reconsidered by Batchelor ,55], among others, in the light of Kolmogorov 1941 
theory 
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Richardson's description of relative dispersion is based on a diffusion equation for the probabiHty density function 
p(r,t) where r{t) = ^2{t) — ^i{t) is the separation of two trajectories generated by (|40p . In the isotropic case the 
diffusion equation can be written as 

dj^l^^l^O. dp(^^ (56) 
ot r'^ or or 

where the turbulent eddy diffusivity was empirically established by Richardson to follow the "four-thirds law" : K{r) — 
koe^/^r'^/^ in which fco is a dimensionless constant. The scale dependence of diffusivity is at the origin of the accelerated 
nature of turbulent dispersion: particle relative velocity grows with the separation. Richardson empirical formula is 
a simple consequence of Kolmogorov scaling in turbulence, as first recognized by Obukhov [56i |. 

The solution of (|56p for (5-distributed initial condition has the well known stretched exponential form 

where A ~ 2187/22407r'^/^ is a normalizing factor. Of course, the assumption the relative dispersion can be described 
by a self-similar process as (fSG]) rules out the possibility of intermittency and therefore the scaling exponents of the 
moments of relative separation 

(r2"(t)} =C2„£"r" (58) 

have the values a„ = 3n, as follows from dimensional analysis. All the dimensionless coefficients C2n are in this case 
given in terms of fco and a single number, such as the so-called Richardson constant C2, is sufficient to parameterize 
turbulent dispersion. 

The hypothesis of self-similarity is reasonable in the presence of a self-affine Eulerian velocity field, such as in the 
case of two-dimensional inverse cascade where the dimensional exponents a2n = 3n/2 have indeed been found [13] ■ 
An analysis of Lagrangian trajectories generated by a kinematic model with synthetic velocity field [5^ has shown 
that Lagrangian self-similarity is broken in the presence of Eulerian intermittency. In this case it is possible to extend 
the dimensional prediction for the scaling exponents by means of the multifractal model of turbulence. 

From the definition of relative separation 

l{rPit)) = {rr'-'5v{r)) (59) 

where dv{r) is the velocity increments between the two trajectories. Using the multifractal representation ([7]) we can 
write 

^(rf(0) J d/irf-i+"+3-^('') (60) 

The time needed for the pair separation to reach the scale r is dominated by the largest time in the process, associated 
to the scale r and therefore given by i ~ r^"''. This leads to 

j^{r'{t)) J (61) 

The integral is evaluated by saddle point method and gives the final result {r^it)) ~ t"'' with scaling exponents 

'p + 3-D{h) 



ar, = inf 

h 



1 



(62) 



From the standard inequality of the multifractal formalism ^ one obtains that even in the presence of intermittency 
a2 = 3. As in the case of single particle dispersion (|50p also here this is a consequence of the presence on the first 
power of £ in (|58|) for n = 1. 

The scaling exponents ap satisfy the inequality ap/p < 3/2 for p > 2. This amounts to say that, as time goes on, 
the right tail of the particle pair separation probability distribution function becomes narrower and narrower. In other 
words, due to the Eulerian intermittency particle pairs are more likely to stay close to each other than to experience 
a large separation. 

The multifractal prediction ([5^ has been checked in synthetic model of fully developed turbulence [s^ where the 
equivalent Reynolds number is very large. In the case of numerical or experimental data, finite Reynolds effects make 
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FIG. 8: Relative dispersion (r^(i)) rescaled with (r^(f))"'"''^ for p = 1,3,4 (+) for Lagrangian simulations in a synthetic velocity 
field corresponding to Rx ~ 10^. The almost constant plateau indicates a relative scaling in agreement with prediction (|62|l . 
For comparison we also plot (r''(t)}_rescaled with the non intermittent prediction {r^(f))^ (x) clearly indicating a deviation 
from normal scaling. Figure from (SS] . 



very difficult to measure the corrections to dimensional exponents. We remark that finite Reynolds effects are more 
important in Lagrangian dispersion than in Eulerian statistics: as a consequence of the accelerate nature of relative 
motion a large fraction of pairs exits the inertial range after a short time. 

To overcome these difficulties in Lagr ang ian statistics, an alternative approach based on exit time statistics has been 
proposed for Lagrangian dispersion [58l.l59j. In close analogy with the exit time approach described in Section FlII C[ one 
computes the doubling times Tp(i?„) for a pair separation to grow from threshold i?„ to the next one Rn+i- Averages 
are then performed over many particle pairs. The outstanding advantage of averaging at fixed scale separation, as 
opposed to averaging at a fixed time, is that crossover effects are removed since all sampled particle pairs belong to 
the same scales. 

Neglecting intermittency, the doubling time analysis can be used for a precise estimation of the Richardson constant 
C2. From the first-passage problem for the Richardson model ([56]) one has 60]: 

^^^^^^^ - 2^^^^^^ 

from which one obtains 

143 - 1)3 i?2 

By using this expression it is possible to estimate from DNS data at moderate Reynolds C2 = 0.50 ± 0.05 [13, [6l| 
which is in agreement with the experimental determination [6^. 

Intermittency effects are evident in higher order statistics of doubling times. In particular, one expects for the 
moments of inverse doubling times, {{l/Tp{R)Y) a power-law behavior 

p 

) - R'^'' (65) 



Tp{.R) 



with exponents /3p connected to the exponents a„ j58| . Negative moments of doubling time are dominated by pairs 
which separate fast; this corresponds to positive moments of relative separation. By using the simple dimensional 
estimate T{R) ^ R/6v{R) one has the prediction 

Pp^Cp-P, (66) 

where are the scaling exponents of the Eulerian structure functions ([5]). 

The multifractal prediction (|66p is found to be consistent with numerical data at moderate Reynolds number. More 
important, as shown in Figure [51 exit-time statistics is sufficiently accurate for discriminating between intermittent 
and dimensional scaling in Lagrangian statistics. 
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FIG. 9: The inverse exit time moments, {[^/Tp{r)Y)^^''^ , for p — 1,...,4 compensated with the Kolmogorov scalings (sohd 
hnes) and the multifractal predictions (dashed hnes). Results from DNS at R\ = 284, initial particle separation ro = 1.2r] and 
ratio between exit times p — 1.25. Figure from [6l| ]. 



VII. DISPERSION IN TWO-DIMENSIONAL CONVECTION: MULTIFRACTAL ANALYSIS OF 

MORE-THAN-SMOOTH SIGNALS 

Thermal convection in two-dimensions provides an example of Bolgiano-Obukhov scaling of turbulent fluctuations. 
Without entering in the details, we recall that within Boussinesq approximation, Bolgiano-Obukhov argument assumes 
a local balance between buoyancy force and inertial term |64| . In the case of two-dimensional turbulence, in the 
presence of a mean temperature gradient, Bolgiano-Obukhov scaling is expected to emerge in the inverse cascade of 
energy with velocity fluctuations given by the scaling law [f35] : 

<5w(r) cx4^^(/3g)2/V3/5, (67) 

where Et is the (constant) flux of temperature fluctuations, (3 is the thermal expansion coefficient and g is the gravity 
acceleration. The prediction (|67p has been checked in both laboratory experiments [66j and in high resolution direct 
numerical simulations [67| which have also shown the absence of intermittency corrections (which is a common feature 
of two-dimensional inverse cascades) . 

We now consider the increments of velocity for Lagrangian tracers transported by Bolgiano turbulence. By extending 
the dimensional argument of Section IIV Al to the general case of velocity scaling exponent h one obtains [62] 

5v{t) oc Vit/T^f/^^-'^'^ (68) 

At variance with Navier-Stokes turbulence, from (|67|) h = 3/5 and therefore q = h/{\ — h) = 3/2 > 1, i.e. velocity 
increments in the inertial range are smoother than Ci signals, the latter denoting the class of differentiable signals. 
This implies that Lagrangian structure functions (j48p are dominated by non-local contributions from the large scale 
L which scale as 

Svit)^TLidtVL){t/TL) (69) 

and therefore give the scaling exponents = p. This set of scaling exponent is trivially universal for any velocity 
field with h > 1/2 and therefore a standard analysis of Lagrangian velocity fluctuations is unable to disentangle the 
non trivial scaling component of the signal [s^l . 

The statistical analysis of more than smooth signals has been recently addressed on the basis of an exit-time statistics 
[68| in which one considers the time increments T{6v) needed for a tracer to observe a change of 6v is its velocity. 
Now, among the two contributions, in the limit of small 6v(t), the differentiable part (|69p will dominate except when 
the derivative dtVL vanishes and the local part (|68|) becomes the leading one. For a signal with 1 < g < 2, its first 
derivative is a one-dimensional sclf-affine signal with Holder exponent ^ — q ~ 1, which thus vanishes on a fractal set 
of dimension D = l — ^ = 2 — q. 

Therefore, the probability to observe the component 0{t'^) is equal to the probability to pick a point on the fractal 
set of dimension D, i.e.: 

P(r ~ Sv^^") ~ T^-^ - {Svy-^/i . (70) 
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By using this probab ility for computing the average p-order moments of exit-time statistics one obtains the following 
bi-fractal prediction [63] 

{TP{5v)) - Sv^'^ , with fip = min(p, ^ + 1 - ^) ■ (71) 

According to prediction (j7ip . low-order moments [p < 1) of the inverse statistics only see the differentiable part of 
the signal, while high-order moments {p > 1) are dominated by the local fluctuations 0{t'^). 

Figure [TUl shows the first moments of exit times {TP{6v)) computed from a direct numerical simulation of two- 
dimensional Boussincsq equation forced by a mean (unstable) temperature gradient which generates an inverse cas- 
cade with Bolgiano-Obukhov scaling. Particles are advected with (|40|) and velocity fluctuations are collected along 
Lagrangian trajectories. The bifractal spectrum predicted by ([71]) is clearly reproduced. We remark that the fact that 
for p > I exponents follows the linear behavior fip — {2p -I- l)/3 indicates the absence of intermittency in Lagrangian 
statistics. This feature is a consequence of the self-similarity of the inverse cascade in two-dimensional Bolgiano 
convection. 
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FIG. 10: The exit time scaling exponents /ip from Lagrangian velocity fluctuations in a DNS of thermal convection with Bolgiano- 
Obukhov scaling. Lines represent the bifractal prediction (|7ip and the error bars on the exponents have been estimated by 
evaluating differences in Hp changing the fitting interval. Figure from [63 |. 



VIII. CONCLUSIONS 

Starting from the seminal work of Kolmogorov we considered the statistical features (mainly scaling properties), 
both Eulcrian and Lagrangian, of the fully developed turbulence in the framework of the multifractal model, i.e. in 
term of D{h). The hard, still unsolved, problem is, of course, to compute D{h) from first principles. Up to now the 
unique doable approach is to use multiplicative models motivated by phenomenological arguments. The non trivial 
result is the fact that, once D{h) is obtained with a fit of the experimental data from the scaling exponents Cp, then 
the other predictions obtained in the multifractal framework, e.g. the pdf of the velocity gradient, the existence of an 
intermediate dissipative range, the scaling of Lagrangian structure functions, are well verified. 
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APPENDIX 



APPENDIX A: SYNTHETIC TURBULENCE: HOW TO GENERATE MULTIAFFINE STOCHASTIC 

PROCESSES 

In this Appendix we describe two methods for the generation of multi-afhne stochastic signals [g^, [t^], whose 
scahng properties are fuUy under controL One is based on a dyadic decomposition of the signal in a wavelet basis 
with a suitable assigned series of stochastic coefficients [6^. The second is based on a multiplication of sequential 
Langevin-processes with a hierarchy of different characteristic times [t^ . 

The first procedure is particularly appealing for modeling of spatial turbulent fiuctuations, because of the natural 
identification between wavelets and eddies in the physical space. The second one looks more appropriate for mimicking 
the turbulent time evolution in a fixed point of the space. 

Using the two methods it is possible to build a rather general {d+ l)-dimcnsional process, w(x, t), with given scaling 
properties in time and in space. 



1. An algorithm based and dyadic decomposition 

A non-sequential algorithm for 1-dimensional multi-affine signal in [0, 1], v{x), can be introduced as [69^1: 

v{x)^Yl — T-^l ^^^^ 

n=l k=l \ f-n / 

where we have a set of reference scales in = 2~" and (p{x) is a wavelet-like function [tH, i.e. of zero mean and rapidly 
decaying in both real space and Fourier-space. The signal v{x) is built in terms of a superposition of fluctuations, 
(fi((x — Xn.k)/in.) of characteristic width and centered in different points of [0, 1], Xn.k = {2k+ l)/2"+^. It has been 
proved [70| that provided the coefficients a„^fc are chosen by a random multiplicative process, i.e. the daughter is 
given in terms of the mother by a random process, an+i,k' — Xo,n,k with X a random number identical, independent 
distributed for any {n, k}, then the result of the superposition is a multi-affine function with given scaling exponents, 
namely: 

{\v{x + e) - v{x)\p) e"^" , 

with ~ — p/2 — log2((ArP)) and £n < ( < In this Appendix ((•)) indicates the average over the probability 
distribution of the multiplicative process. 

Besides the rigorous proof, the rationale for the previous result is simply that due to the hierarchical organization 
of the fluctuations, one may easily see that the term dominating the expression of a velocity fluctuation at scale R, 
in (|A1[) is given by the couple of indices {n, k} such that n ~ log2{R) and x ~ Xn,k, i-e. v{x + i) — v{x) ^ Un^k- The 
generalization (|A1[) to d-dimension is given by: 

n=l fc=l \ / 

where now the coefficients {an,k} are given in terms of a d-dimensional dyadic multiplicative process. 



2. A sequential algorithm 

Sequential algorithms look more suitable for mimicking temporal fluctuations. With the application to time- 
fluctuations in mind, we will denote now the stochastic 1-dimensional functions with u{t). The signal u{t) is obtained 
by a superposition of functions with different characteristic times, representing eddies of various sizes JO] : 

N 

u{t) = J2un{t) . (A2) 

Ji=l 

The functions Un{t) are defined by the multiplicative process 



Un{t) = gn{t)xi{t)X2{t) ■ ■ ■ Xn{t) , 



(A3) 
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where the gn(t) are independent stationary random processes, whose correlation times are supposed to be r„ = (in.)", 
where a — 1 — h (i.e. t„ are the eddy-turn-over time at scale in) in the quasi-Lagrangian frame of reference [72| and 
a = 1 if one considers u{t) as the time signal in a given point, and (g^) = (in)'^^, where h is the Holder exponent. For 
a signal mimicking a turbulent flow, ignoring intermittency, we would have h = 1/3. Scaling will appear for all time 
delays larger than the UV cutoff rjv and smaller than the IR cutoff ri. The Xj(t) are independent, positive defined, 
identical distributed random processes whose time correlation decays with the characteristic time tj . The probability 
distribution of Xj determines the intermittency of the process. 

The origin of (|A3|) is fairly clear in the context of fully developed turbulence. Indeed we can identify m„ with the 
velocity difference at scale in and Xj with (ej/ej-i)^^^, where Sj is the energy dissipation at scale £j |70 |. 

The following arguments show, that the process defined according to (|A2IA3[) is multi-afSne. Because of the fast 
decrease of the correlation times tj — (ij)"', the characteristic time of is of the order of the shortest one, i.e., 

Tn = {in)"- Therefore, the leading contribution to the structure function Sq{T) = {\u{t + r) — u{t)\'^) with t ^ Tn 
stems from the n-th term in (jA2|) . This can be understood noting that in u(t -\- t) — u(t) — 'Y^^=i['^k{t + t) — Uk{t)\ 
the terms with k < n are negligible because Uk{t + t) ~ Uk{t) and the terms with k > n are sub-leading. Thus one 
has: 

hq log2((^''>) 

S,{rn)^{\Unn^{{\gnn){{x'^)r^T;F (A4) 



and therefore for the scaling exponents: 



a a 



The limit of an affinc function can be obtained when all the Xj are equal to 1. A proper proof of these result can be 
found in [t^. 

Let us notice at this stage that the previous "temporal" signal for a — 1 — /i is a good candidate for a velocity 
measurements in a Lagrangian, co- moving frame of reference (7^ . Indeed, in such a reference frame the temporal 
decorrelation properties at scale are given by the eddy-turn-over times r„ = {in)^~^ ■ On the other hand, in the 
laboratory reference frame the sweeping dominates the time evolution in a fixed point of the space and we must use 
as characteristic times of the processes a;„(t) the sweeping times rl"' = £„, i.e., a = 1. 
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